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I.  INTRODUCTION  AND  PRELIMINARIES 


This  report  is  intended  as  the  first  of  a  series  in  which  the  problem 
of  heat  transfer  at  a  solid  boundary  will  be  mathematically  investigated 
with  the  ultimate  objective  of  developing  accurate  and  efficient  numeri¬ 
cal  algorithms  to  solve  the  appropriate  governing  equations.  The  basic 
line  of  attack  will  be  to  use  the  tools  of  asymptotic  analysis  to  obtain 
approximate  solutions  for  equations  in  the  form 

ut  =  aO)  uxx  +  b(x)  +  c(x)  u  (1.1) 

valid  for  small  times  (in  a  non-dimensionalized  sense  to  be  made  more 
precise) .  In  particular  these  will  be  obtained  through  the  use  of  a 
"Diffusion  Equation  Solution  Sequence"  or  DESS,  a  concept  which  was 
introduced  in1.  The  approximate  solutions  will  then  be  incorporated 
into  a  numerical  scheme  which  advances  over  longer  time  intervals  by 
repeating  the  basic  short-time  algorithm  as  frequently  as  necessary. 

A  number  of  problems  of  engineering  concern  can  be  treated  using  this 
technique.  The  problem  of  particular  importance  is  that  of  heat  transfer 
in  gun  barrels.  More  generally  any  "transport"  equation  such  as  discussed 
in2  can  be  treated  in  similar  fashion. 

In  the  present  report  we  shall  not  accomplish  these  longer  term 
objectives  but  merely  wish  to  lay  the  foundations  by  defining  and  investi¬ 
gating  the  properties  of  certain  special  functions  to  be  denoted 

H*(x,t)  and  vn(x,t) .  These  solve  the  diffusion  equation  and 

satisfy  certain  special  initial  conditions.  In  Section  VI  we  shall  show 
how  these  functions  can  be  effectively  utilized  for  obtaining  asymptotic 
expansions  for  solutions  of  initial-boundary  value  problems  for  the 
diffusion  equation,  valid  for  small  times. 


The  number  of  independent  variables  in  diffusion  problems  can  frequent¬ 
ly  be  reduced  from  two  to  one  by  considering  the  similarity  parameter 
-1/2 

xt  ,  t>0„  This  is  true  in  the  present  analysis  and  we  have  found  it 

convenient  not  to  introduce  the  functions  H*  H  and  v  directly  but 

Y  Y  n  7 


first  to  consider  the  related  functions  H  ,  H*  and  v 

Y  Y  n 


depending  only 


-1/2 

on  a  single  variable  z  (which  may  be  taken  as  x  t  ) .  The  first  two 
functions  are  generalizations  of  the  well  known  repeated  integrals  of 
the  error  function  and  the  last  is  a  related  polynomial. 


1J'.  F.  Polk3  "Asymptotic  Expansions  for  the  Solutions  of  Parabolic 
Differential  Equations  with  a  small  parameter "3  Ph.D.  dissertation 3 
Department  of  Mathematics 3  University  of  Delaware 3  Newark 3  DE,  1977. 

2P.  J .  Roache 3  "Computational  Fluid  Dynamics "3  paqe  18,  2nd  Edition* 
Hermosa  Publishers ,  Albuquerque 3  NM3  1976. 
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Except  for  Section  VI  the  contents  of  this  report  represent  a  simpli¬ 
fied  version  of  Chapter  1  of  the  author's  doctoral  dissertation1.  Some 
of  the  results  obtained  therein  have  been  re-presented  here  in  report 
form  mainly  to  provide  a  conveniently  available  reference  for  subsequent 
work. 


The  following  notational  conventions  will  be  used  in  this  report 

f*  (z)  =  f  C-z) 

f*  (x,t)  =  f  (-x,t) 

f'  -  h  f 

fCn)(z)  ,  ^  f  (!) 

dz 

Dn  f(x,t)  =  \  f  (x,t) 

9x 

D  f(x,t)  =  —  f  (x,t) 

at 


[a]  =  greatest  integer  a 

R  =  (-00,  °°)  =  all  real  numbers. 


Because  the  gamma  function  will  be  frequently  encountered  in  the 
discussion  it  is  convenient  to  use  the  more  simplified  notation  of  the 
factorial  function.  This  is  defined  by 

oo 


for  a  >-l  and  by 


a! 


(a+1) ' 
a+1 


for  a  <-l,  a  ^  -2,  -3,.... 


This  relates  to  the  usual  gamma  function  by 


al  =  r  (a+1) 


The  reciprocal  of  the  gamma  function  is  known  to  be  entire  with  zeros  at 
0,  -1,  -2,  ...;  thus  the  function  1/a!  is  well  defined  and  finite  for 
all  aeR  with 


1_ 

a.' 


0 


(1.2) 
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for  a  =  -1,  -2,  ...  .  Using  factorial  notation  the  well  known  duplica¬ 
tion  formula  for  the  gamma  function3. becomes 


2a  /r 

a!  "  ( (a-l)/2) !  (a/ 2)1  (1.3) 

II.  THE  FUNCTIONS  H 

Y 

For  any  YeR  define  the  "canonical"  jump  functions 


0 

x<0 

hY  w  = 

xY/y'. 

x>0 

(2.1) 

and 

h* O)  = 

hY  C'x)i 

> 

for 

Y  =  -1 

define 

“-1 

(z)  =  (4tt)  '*'//2  exp 

(-z2/ 4) 

(2.2) 

and 

for  y> 

-1  let 

H 

Y 

(z)  =  fh 

R 

,  (s)  H_x  (s- 

■z)  ds 

(2.3) 

where  R  =  (-00,  °°) ;  in  particular  note  that 


H0  (z)  =  (4tt)"1//2  J  exp  (- (s-z) 2/ 4)  (2.4) 

o 

=  erfc  (-z/2)  . 

More  generally,  comparing  with  equation  (7.2.3)  of3’  page  299. 

00 

Hn  O)  =  (4tt)  1//2  / s”  exp  (-  (s-z) 2/ 4)  ds 

o  m!  (2.5) 

„n-l  .n  _  , 

=2  i  erfc  (-z/2) 


%.  Abramowitz  and  I.  A.  Stegun,  editors ,  "Handbook  of  Mathematical 
Functions  with  Formulas ,  Graphs  and  Mathematical  Tables ",  p  256,  National 
Bureau  of  Standards  Applied  Mathematics  Series  #55,  Second  Edition,  1964. 
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Where  in  erfc  denotes  the  n-  th  integral  of  the  complementary  error 


function.  For  y  <  -1  w.e  define  H^,  by 

_  _(n) 

H  (z)  =  H  (z) 

Y  Y+n  v 


(2.6) 


where  the  superscript  denotes  differentiation  n  times  and  n  =  [-y]  = 
integer  part  of  -y.  In  or_der  that  this  definition  for  y  <  -1  be  meaning¬ 
ful  it  is.  necessary  that  H  be  differentiable  for  -1  <_  y  <  0*  This 
is  a  consequence  of  the  following  proposition. 

Proposition  1;  The  functions  H  are  entire  functions  for  all  y  e  R  and 

J»  _ 

H  =  H  (2.7) 

Y  Y-n 

for  any  non-negative  integer  n. 

Proof:  The  function  H  ^  is  clearly  entire.  For  y  >  -1  consider  the 
related  functions 


(z)  =  (4tt) ly/2  y'.  exp  [z2/4]  (z) 


■fi 


-J  exp  [  -  (s  -  2sz)/4]  ds 
o 

These  functions  are  clearly  positive  and  have  n  -  th  derivatives  given  by 

oo 

=  2  n  J  s"+n  exp  [  -  (s2  -  2sz)/4]  ds 

(z) 


0-n 

=  2  g 


Y+n 


>  0. 


Thus  gy(z)  is  absolutely  monotonic  on  the  real  line.  Such  functions 
are  necessarily  analytic4.  We  may  conclude  that 

Hy  (z)  =  (4ir)  “V2  exp  (-z2/4)  gy(z) 

is  entire  for  y  >  -1  also.  To  show  that  (2.7)  holds  for  n=l  and  y  >  0 
we  may  differentiate  (2.3)  to  obtain 


H"  (z)  =  H  (z) 

Y  Y_i 


^D.  V.  Widder,  "Tne  Laplace  Transform  ",  page  146  3  Princeton  University 
Press,  Princeton,  NJ  1941. 
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But  from  (2.2)  and  2.4)  we  see  that  this  must  hold  for  y  =  0  also  and 
can  be  extended  to  all  7  using  definition  (2.6).  For  n  >_  2  equation 
(2.7)  then  follows  by  simple  induction.  This  concludes  the  proof. 

The  present  version  of  the  argument  showing  that  H  is  entire  is 
due  to  Mr.  W.  0.  Egerland  of  the  Ballistics  Research  Laboratory5. 

A  useful  recursive  formula  can  now  be  obtained.  For  y  >  0  we  may 
differentiate  (2.3)  under  the  integral  sign  and  obtain 


H  i  (z) 
y-1  v 


h;  (z) 

I  As  -  ■)  hY  (s)  (s  -  z)  ds 

R 

hTLf K+i  (s)  h_i  (s'z)  ds 

R 

-  j  /hy(s)  H_1  (s-z)  ds 
R 


Lull 

2 


Hy+1  (z) 


(z) 


By  shifting  the  index  and  rearranging  terms  the  last  equation  can  also 
be  written  in  either  of  the  forms 


Y  (z) 


H 

Y 


(z) 


z  (z)  +  2  Hy_2  (z) 

(2.8a) 

LlLLL  h  (z-)  z  jp  (t-) 

2  y+2  lJ  2  y+1  LZJ * 

(2.8b) 

By  differentiating  these  formulas  it  is  easy  to  extend  their  validity 
to  all  choices  of  y. 

Let  us  now  derive  a  power  series  expansion  for  (z) .  First 

note  that  for  y  >  -1,  using  standard  tables  of  definite  integrals6 
we  have 

5&L  0.  Egerland ,  Private  Communication ,  US  Army  Ballistic  Research  Labora¬ 
tory,  Aberdeen  Proving  Grounds,  MO,  1977. 

6  -  j  Handbook  of  Mathematical  Tables,  Supplement  to  Handbook  of  Chemis¬ 

try  and  Physics,  page  324,  Second  Edition,  The  Chemical  Rubber  Company, 
Cleveland,  OH,  1964. 
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Hy  (0)  =  (4tt]) -1/2  J  1_  exp  (-s2/4)  ds 

0 

-  2lif:>—  /r cy-d/2  exp  (.r)  dr 
Y  •  o 

oY  r4^-l/2 

=  2_I4rrJ -  ((y-l)/2)I 

T  • 

But  then  from  the  duplication  formula  (1.3) 


(2.9) 


H  (0) 
Y 


1 

2  (y/  2) 


(2.10) 


This  can  be  shown  to  be  valid  for  y  <  -1  also, through  the  use  of  recursion 
formula  (2.8).  It  follows  that 


CO 

H  (z)  =  £  H  (n)(0)  zn  /  n! 

Y  n=o  Y 


=  V  H  (0)  zn  /  n' 

Z-*  Y-n 


n=o 


Y 


E 


n 


2  ^  ( (Y-n)/ 2) !  nl 
n=o 


which  converges  for  all  z  since  is  entire. 


(2.11) 


A  different  expansion  for  H^(z)  is  useful  for  z  <  0  since  it  turns  out 

that  H  (z)  vanishes  exponentially  as  z  +  -  Recall  from  the  proof  of 
Proposition  1  that  for  y  >  -1 


Hy  (z)  =  H_!  Cz)  gy  (z)  /Y! 

Because  the  function  g  (z)  is  entire  it  may  be  expressed  as  a  power 
series.  Thus  ^ 


gY  =  Z)  gy^CO)  zn  /  n! 
n=o 
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00 


which,  for  (2.9), 


consequently 


H 

Y 


(z) 


=  E  2"n  gy+n  (0)  zn  /  n! 
n=0 

=  C4ir)-l/2£  2_n  (Y+n)  I  H  (0)  zn  /n! 
n=o  1 


=  2Y  V  ((Y+n-l)/2) !  Zn 

^  n! 


n=o 


H ,  Cz)  2  zn 

n=o  1  *  n< 


(2.12) 


Using  induction  this  formula  can  be  extended  to  all  y  except  y  =  -1, 

-2,  -3,  ...  where  the  expression  for  the  coefficients  becomes  indeter¬ 
minate  and  requires  special  interpretation.  We  shall  obtain  an  alterna¬ 
tive  formula  for  these  cases  presently  (equation  (2.20)  below). 


To  simplify  the  evaluation  of  the  functions  H  note  that  once  H  (z) 
(z)  have  been  determined  for  some  y  then^the  recursive  ^ 

formula  (2.8)  can  be  used  to  evaluate  H  n  (z)  for  any  integer  n.  In 

particular  we  can  write 


Vn  «  ■  PY,n  O)  Hy  Cz)  *  qy>n  (z)  (z)  (2.13) 

where  p^  ^  (z)  and  q^  ^  (z)  are  polynomials -of  degree  [n+1 1  and  |n|-l 
respectively  which  satisfy  the  following  recursion  formulas 


P  (z)  =  1 

rY,o  v  J 

Yo  (z)  =  0 

P  t  (z)  =  0 
*Y>-1  ^ 

Y-i  (z)  =  1 

(Y+n)  P^n(z)  = 

z  P  .  (z)  +  2  p  _ 

y,n-l  v  ry,n-2 

(z) 

(2.14) 

(Y+n)  q^n  (z)  = 

z  q  .  (z)  +  2  q 

ny,n-l  J  ny,n-2 

(z) 

(2.15) 

These  may  be  verified  without  difficulty  using  induction  and  we  will  not 
include  the  details  here.  The  first  few  polynomials  are 
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’y,-4  (z)  =  i  Y  z2  +  T  Y  (Y"2) 

qy,_4  =  -  i  z3  -  j  c2^"3)  z 

y,-3  M  ’  -TY  2 

f  -v  1  2  1.  , . 

V-3  (Z)  =  4  2  +  2  CY_1) 

Y>  “2  ^  =  2  Y 

qY,-2  (Z)  =  _  I  2 

o 

II 

v _ / 

r— 1 

1 

•v 

>- 

qY,"l  (Z)  =  1 

'y.o  w  -  1 

qY,0  (Z)  =  °  (2 

Y»1  (Z)  '  Y*1 

qY,l  (Z)  (Y+l) 

,  f2l  .  k  *  ihiil 

Y,2  '•Zj  (y+1)  (Y+2) 

r  2z 

y >2  (-Zj  (y+1)  (Y+2) 

PY,3  M 


z+2  (2y+3)  z 

(Y+l)  (y+2)  (y+3) 


q  (z)  =  ?-z  +  4  CY+2J 

4Y,3  (Y+l)  (Y+2)  (y+3) 


Of  particular  interest  are  the  functions  H  (z)  where  n  is  an 


integer.  Choosing  y=0  in  (2.13)  yields. 


n 


H  (z)  =  p  (z)  H  +  q  (z)  H  _ 

n  o,n  o  no,n  ^  '  -1 

where  from  (2.2)  and  (2.4) 

hq  ■  j  erfc  (-z/2) 

and 

H  *  (4ir)  ^2  exp  (-z2/4) 

Thus  from  (2.16) 


H-4 

(z) 

=  (-1/2) 3  (z3  -  6z)  H_1 

ff-3 

(z) 

=  (1-1/2) 2  (z2  -  2)  H_1 

H-2 

(z) 

=  (-1/2)  z  H_1 

(2.17) 
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(2.18) 


\  (z)  = 

h2  (z)  -  [ 


H3 


z  H  +  2  H  . 
o  -1 


(z  +  2)  H  +  2  Z 
v  J  o 


/  2 


■[ 


(z3  +  6z)  Hq  +  2  (z2+4)  " 


H.i]  /S! 


h4  Cz) 


[(z4  +  12z2  +  12)  Hq  +  2  (Z3  +  10z)  H_1J  /4' 
The  cases  n=  -1,  -2  simplify  since  it  can  be  shown  using  (2.14) that 


P  (z)  =  0 

o,n  v 


(2.19) 


for  n  <  -1  and  thus  from  (2.17) 


H 


n  (z) 


=  „  (z)  (z) 


o,n 


(2.20) 


n=  -1.  -2.  ...  .  For  Ini  <  12  the  polynomials  p  and  q  are  listed 
’  ii_  L  o,n  o,n 

explicitly  in  the  Appendix. 

The  asymptotic  behavior  of  H  (z)  as  |  z  |  -*  00  is  completely_determined 
for  n=  -1,  -2,  by  the  representation  (2.20).  The  behavior  of  (z) 

for  y  =  ~1>  -2,  ...  is  somewhat  different,  especially  for  z  ■*  +  °°. 

We  conclude  this  section  by  stating  without  prooif1  two  propositions 
which  characterize  the  asymptotic  properties  of  H^. 

Proposition  2:  For  any  interger  N  >_  0  and  any  y  eR  with  y=  -1,  -2,  ... 

there  exists  a  constant  K  XT  >  0  such  that 

y,N  - 


H  (z)  =  (-2/z)Y+1  H_1  (z) 


LW  t-2)"  *  w 


n=o 


(2.21) 


where  R  „  fz)  is  a  remainder  term  satisfying 

y,N  J 


R  (z)  <  K  z| 

Y»N  ^  J  ~  Y >N  1  1 


■2(N+1) 


uniformly  for  z  <  0  . 
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Proposition  3;  For  any  integer  n  ^  0  ,  y  e  R  and  zQ  >  0  there  exists  a 


(2.22) 


uniformly  for  z  ^  z^. 


Note  that  Proposition  3  remains  valid  when  y  =  -1,-2,...  although 
it  is  somewhat  degenerate  in  these  cases  since  the  summation  term 
vanishes  in  view  of  (1.2).  Proposition  2  can  also  be  shown  to  hold 
when  y  is  a  negative  integer  by  properly  interpreting  the  indeterminate 
expression 


(y  +  2n) ! 

Y* 


III.  THE  FUNCTIONS  H*  and  v 

Y  i 


n 


From  the  notation  introduced  in  Section  I  we  have 


H*  (z)  =  H  C-z) 

Y  Y 


(3.1) 


for  all  y  e  R.  From  (2.7)  and  (2.8)  it  follows  immediately  that  for 
any  y  e  R 


H*  (z) 
Y 


(3.2) 


YH*  (z) 


-z  H*  .  (z)  +  2  H*  _  (z) 

y-1  J  y-2  *■  J 


(3.3a) 


and 
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(3.3b) 


(z)  =  ( (y+2) / 2)  H*  (z)  +  £H*  (z) 

Y  Y*2  2  Y+1 

From  (2.2)  it  is  clear  that  H  =  H*  and  thus  from  (3.2)  and  (2.7) 
we  have  ~ 

H*  (z)  =  (-l)n+1  Hn  (z)  (3.4) 

for  n  =  -1,-2, -3, . . . 

The  functions  are  defined  for  integer  values  of  n  by 


v»  (z)  '  /  ht  H-i  ts-z)  dz  (3.5) 

Comparing  with  (2.3)  we  obtain  the  identity 

7n  (z)  =  ffn  (z)  +  C-Dn  H*  (z)  (3.6) 

which,  from  (3.4),  gives 

vn  (z)  =  0  (3.7) 

for  n  »  -1,-2, -3, . . . ;  from  (2.7)  and  (3.2)  we  see  that 

Vn(k‘)  (z)  =  7_k  (z)  (3.8) 

and  from  (2.8a)  and  (3.3a)  we  have 

n  7n  (z)  =  z  7n_1  (z)  +  2  vn_2  (z)  (3.9) 

Evaluating  the  integral  in  (3.5)  explicitly  yields 

\  (z)  =  1  (3.10) 

V1  (z)  =  z 


But  then,  comparing_(3 .9)  with  (2.14)  and  (3.10)  with  (2.16)  we  see 
that  the  functions  v^  coincide  with  the  polynomials  p^  ^ 

\  (z)  ■  Po,n  <*)  0.11) 
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for  3ny  n.  Moreover  the  recursion  relation  (3.9)  can  be  used  to  verify 
the  following  representation  for  v^ 


v 

n 


[n/2] 


k=o 


n-2k 
z _ 

k!  (n-2k) ! 


(3.12) 


n  =  0,  1,  2,  ...  where  [n/2]  =  greatest  integer  £  n/2.  For  n  _<  12 
these  are  listed  explicitly  in  the  Appendix.  Equation  (3.10)  can 
also  be  inverted  to  obtain 


[n/2] 

=  ^2  (-Dk  V2k  (Z) 

i  __  k! 


(3.13) 


IV.  THE  FUNCTIONS  H 

Y 

In  this  section  the  relationship  between  the  functions  discussed 

in  Section  II  and  the  diffusion  equation  is  made  clear  by  introducing 
the  functions  H^.  These  are  defined  for  any  y  e  R  by 

Hy  (x,t)  =  yfty  Hy  (x/Vt  )  (4.1) 

for  t  >  0,  with 

Hy  (x,0)  =  h^  (x)  (4.2) 

for  t  =  0.  Specifically  for  y  =  -1  and  y  =  0  we  have 

H_x  (x,t)  =  (47rt) - ^  exp  (-x2/4t)  (4.3) 

HQ  (x,t)  =  j  erfc  (- x/2\/t )  (4.4) 

For  y  >  -1  these  functions  can  also  be  given  an  integral  representation 
using  (2.3) 


Hy  (x,t) 


(s)  H 


(s-x,t) 


o 


ds 


[-(s-x)2/4t]  ds 


(4.5) 
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The  assignment  of  initial  values  in  (4.2)  is  not  arbitrary  but 

provides  a  continuous  extension  of  from  t  >  0  into  t  >  0 ,  except 

possibly  at  the  point  x  =  0,  t  =  0.  To  see  this  note  from  (2.21)  that 
for  x  <  0 


iim  H  (x,t)  =  lim  \/tY  H  (x/Vt) 

t+o  '  t->o  Y 

=  lim  \T :Y  exp  (-x2/ 4t)  (-V4t/x)Y+1  I"  I—  +  0  (t/x2)l 
t->o  Lm!  •  J 

=  0 

and  from  (2.22)  for  x  >  0 

iim  H  (x,t)  =  lim  ^ — C.V [1  +  0 (t/x2) ] 
t+o  y  t->o  y ! 


Yl  ' 

The  functions  H  can  be  shown  to  be  continuous  at  (0,0)  if  and 
only  if  y>0  and  bounded  if  and  only  if  y  >_  0 .  When  y  is  a  negative 
integer  the  initial  values  (4.2)  are  seen  to  vanish  because  of  (1.2). 
This  does  not  fully  convey  the  limiting  behavior  of  H  as  t->o  since 
this  can  be  properly  expressed  only  in  terms  of  distributions  or 
generalized  functions.  A  more  accurate  formulation  would  be 

H_n  (*,0)  =  6Cn_1)  (x) 

for  n  -  1,2,3,  ...  where  6  (x)  is  the  usual  Dirac  delta  function  and 

<5  indicates  its  k-th  generalized  derivative.  Since  these  concepts 
will  not  be  required  in  our  analysis  we  shall  not  discuss  them  further. 

Evaluation  of  for  the  important  cases  where  y  is  an  integer  can 
be  accomplished  using  (2.17)  from  which  we  have 


Hn  (*,t) 


=  \TtU 


P0,nCx/^  H0(X>V  +V*  Vn  Cx/\fa  H_iO,t) 


(4.6) 


where  the  polynomials  p  and  q  are  given  in  the  Appendix.  In 

u  ^  ii  o  *  rw 

particular  the  first  few  of  these  are 


H-5  =  [(x4  _  12x2  z  +  120  /  C-2t)4]  H_1 
H_4  =  [(x3  -  6xt)  /  (-2t)3]  H_1 
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(4.7) 


H 


H 


-2 


[(x2  -2t)  /  (-2t) 2 J  H  x 
-(x/2t)  H 


H, 


H, 


H 


H 


-1 

x  H  +  2t  H 


-1 


(x2  +  2t)  Hq  +  2xt  H_1  J  /2 

(x3  +  6xt)  Hq  +  2  (x2t  +  4t2)  H  i  J  /3* 

(x4  +  12x2t  +  12t2)  Hq  +  2  Cx3t  +  10xt2)  H_1J  /4l 


The  properties  of  the  function  follow  directly  from  those  of 
the  functions  H^0  Setting  x  =  0  in  (4.1)  and  using  (2.10)  we  see 


yfs 

V0^  =  ~1WWT 


=  2  \/2  1  -  °‘ 


(4.8) 


Differentiating  (4.1)  and  recalling  (2.7)  and  (2.8)  yields  the 
formulas 


Dt 

H 

Y 

t  >  0 

(4.9) 

D 

X 

H 

Y 

■  Vi 

t  >  0 

(4.10) 

Dn 

X 

H 

Y 

=  H 

Y-n 

t  >  0 

(4.11) 

and 


It  immediately  follows  that  is  a  solution  to  the  diffusion  equation 

[0t  -  D2]  Hy  (x,t)  =  0  (4.12) 

in  the  region  t  >  0.  More  generally  for  any  a  >  0  and  any  xq  the 
function  (x-xq,  at)  satisfies  the  equation 


[d  -  aD2]  H 
L  t  xj 


y  (x"x0>  at)  =  0 


for  t  >  0  and  has  the  initial  values 


Hy  (x-V0)  =  \  (x~V 

From  the  recursion  formulas  (2.12)  we  obtain 


YH  =*H  .  +  2t  H  „ 

Y  y~1  y_2 


H 


Cl/2t)  [  CY+2)  H7+2  -XHytl] 


(4.13) 

(4.14) 

(4.15a) 

(4.15b) 
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for  any  y  e  R. 


A  number  of  bounds  on  the  growth  of  the  functions  H  have  been 

•  1  Y 

established  in1  but  for  present  purposes  we  need  only  the  following: 

For  any  xq  <  0  and  any  T  >  0  there  exists  a  constant  K  >_  q  such  that 

Ihy  (x>t) I  1  K  exp  (-x2/ 4t)  (4.16) 

uniformly  for  all  x  <  x  and  0  <  t  <  T.  This  can  be  verified  by 

“  o  _  - 

noting  from  (2.21)  that  as  y?/f >  00 
Hy  (x,t)  =^tY  Hy  (x/y/t) 

=  (4tt)  (-2^t/x)Y+1  exp  (-x2/4t)  1+  0(t/x2)J 

and  -x  /V  t  i  -x0/V"t  > 0  in  the  region  of  interest. 


V. 

THE  FUNCTIONS 

H*  and  v_ 
Y  n 

The  functions  H*  are 
Y 

defined  for  y  e 

R  by 

H* 

Y 

(x,t) 

=  Hy  (~x,t). 

t  > 

0 

(5.1) 

From  the  properties  of  we  immediately  have 

H* 

Y 

(x,0) 

=  h*  (x) 

(5.2) 

H* 

Y 

(0,t) 

2  hy/2  ^ 

t  >_ 

0 

(5.3) 

Dt 

H* 

Y 

=  h;-2 

t  > 

0 

(5.4) 

D 

X 

H* 

Y 

=  -H*  , 

Y-l 

t  > 

0 

(5.5) 

Dn 

X 

H* 

Y 

=  (-l)n  H* 

Y-n 

(5.6) 

K 

-  aD2 

)  (X'V  at) 

=  0  a>0 

,xeR,t>0 

(5.71 

yH* 

Y 

=-XH 

*  .  +  2t  H*  , 

Y-l  Y-2 

t  >0 

(5.8a) 

H* 

Y 

-  (l/2t)  [  (-,.2)  H*+2  «hO 

t  >  o. 

(5.8b) 
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From  (3.4)  we  have 


H*  =  (-l)n+1  H 

n  ^  J  n 

for  n  =  -1,  -2,  . . . 

The  functions  v  are  defined  for  integer  values  of  n  by 


(5.9) 


n 


v  (x 
n  v 


,t)  =  yftn  vn  (x/Vt) 


for  t  >  0,  with 


vn  (x,Q)  =  xn/n.' 


(5.10) 


(5.11) 


for  t  =  0.  Using  (3.5),  (3.6)  and  (3.12)  we  see  that  these  may  also 
be  written  in  the  forms 


and 


vn  (x,t) 


v  (x,t) 
n  v  ' 


\  (x,t) 


=  /4  H 

R  n'  - 1 


(x-s,t)  ds 


.  n 


Hn  (x,t)  +  (-1)“  H*  (x,t) 


[n/2] 

E  x 


k=o 


n-2k  t 
(iv-2k)  !  kT 


(5.12) 

(5.13) 

(5.14) 


For  n  >_  0  these  functions  are  thus  identical  (except  for  a  factor  nl) 
with  the  well  known  heat  polynomials  discussed  in  Ref  7  and  8;  they  are 
fundamentally  important  because  they  are  the  polynomial  solutions  of 
the  diffusion  equation 


D  -  D 
t  x 


v  (x,t)  =  0 
n  v  J 


(5.15) 


with  the  polynomial  initial  values  (5.11).  Evaluation  along  x  =  0 
clearly  yields 

n/2 


vn  (0,t)  = 


(n/2) ! 
0 


if  n  =  0,2,4,  ... 
if  n  =  1,3,5,  ...  . 


(5.16) 


For  n  <_  -1  note  from  (3.7)  that 


7P.  C.  Rosenbloom  and  D.  V.  W-idder ,  "Expansions  in  Terms  of  Heat 
Polynomials  and  Related  Functions , "  Transactions  of  the  American 
Mathematical  Society ,  92  (1959) ,  pp.  220-266. 

8D.  V.  Widder3  "The  Heat  Equation3  "  Academic  Press ,  Hew  York,  1975. 
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(5.17) 


vn  (x,t)  =  0 

identically,  for  all  x  and  t. 


VI.  INITIAL-BOUNDARY  VALUE  PROBLEM 

As  indicated  in  the  introduction  we  shall  eventually  use  the  functions 
Hy  and  vn  to  develop  asymptotic  expansions,  accurate  for  small 

times,  for  solutions  to  equations  such  as  (1.1).  This  will  be  accom¬ 
plished  in  later  reports.  We  are  already  in  a  position  to  do  this  for 
the  diffusion  equation  itself,  however,  and  shall  now  show  how  various 
boundary  and  initial  conditions  can  be  effectively  handled. 


For  illustrative  purposes  let  us  first  consider  two  simple  but 
important  problems  for  the  semi-infinite  domain  x  >  0,  t  >  0 


[°t  -  aD^]  u 

=  0 

x,  t  >  0 

U  (x,0)  = 

h  (x-x  ) 

Y  o' 

x  >  0 

(BVP) 

u  (0,t)  = 

0 

t  >  o 

and 

[Dt  -  aDx]  " 

=  0 

x,t  >  0 

(BVP) 

w  (x,0) 

=  n 

x  >_  0 

w  (Q,t) 

=  h  (t) 

t  >  0  . 

These  problems 

have  the  exact  solutions 

u  (x, t)  =  H  (x-x 

Y  o 

,  at)  -  H* 
Y 

(x+xq,  at) 

(6.1) 

and 

w  (x,t)  =  2a~Y 

H*y  (x,at) 

(6.2) 

respectively,  as  can  be  verified  by  noting  (4.13)  and  (5.7)  and  by 
directly  substituting  (4.14),  (5.2)  and  (5.3)  into  (6.1)  and  (6.2). 

We  can  also  guarantee  that  these  solutions  are  unique  by  imposing 
additional  growth  conditions  on  u  and  w  but  this  is  not  crucial  to 
our  discussion.  One  important  point  to  note  regarding  these  functions 
is  that  they  are  defined  not  just  for  x  >_  0  but  for  all  x.  Thus, 
they  can  also  be  considered  as  solutions  to  the  following  Cauchy 
or  pure  initial  value  problems  respectively: 
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0 


and 


u  (x,0) 


h  (x-x  ) -  h* 
y  *■  o  y 


(x+xq) 


(IVP)X 


w  (x,o)  =  2a"Y  h*Y  (x)  (IVP)2 

This  equivalence  between  initial-boundary  value  problems  and  Cauchy 
problems  is  a  basic  property  of  the  diffusion  equation;  it  can  also  be 
used  for  dealing  with  the  more  complicated  problems  involving  a  bounded 
domain  0  <  x  <  l  as  we  now  intend  to  show. 


Consider  the  general  initial-boundary  value  problem 


aD  u  =  0  0<t<T,  0<x<& 


u  (x,0)  =  f(x)  0  <  x  <  ^ 


(6.3) 

(6.4) 


u  (0,t)  =  g  (t)  0  <  t  <  T 


(6.5) 


u  (A,t)  =  g  (t)  0  <  t  <  T  (6.6) 

It  may  be  supposed  that  the  functions  f(x),  go(„t)  and  g^(t)  arise  from 

thermocouple  measurements  and  do  not  have  any  specific  "analytic" 
form.  However,  by  using  a  combination  of  polynomials  and  jump  functions 
h  ,  they  can  usually  be  effectively  approximated  with  relatively  few 

terms.  For  instance  we  may  have  the  initial  and  boundary  conditions 


and 


*(x) 


nr  n 

E[v '•£ 

k=o  L  j=o 


bjkhk 


(x 


'V 


e(x) 


gQ  (t)  =  0  =  g1  (t) 

where  e(x)  is  an  error  term  combining  the  errors  in  measurement  and  the 
errors  in  representation  which  is  bounded  in  the  form 

|  e  (x)  |  _<  e^  0  <  x  <  2 

for  some  "small"  number  e  >0.  Using  linearity  the  solution  to 
problem  (6.3)  -  (6.6)  can°then  be  written  in  the  form 
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n  n 

u  cx>t)  =2kuk  <*•*>  *Ev  »,k  tx>t) 

k=o  L  j=o 


+  E  (x,t) 


where  and  denote  solutions  to  (6.3)  which  vanish  along  x  =  0 
and  x  =  H  and  satisfy 

uk  (x,0)  =  xk  Ox 

wjk  (x»°)  =  hk  0-xj) 

respectively.  The  term  E  is  an  error  term  which  will  also  satisfy  the 
bound 

|E  (x,t) |  £  eQ 

because  of  the  maximum  principle  for  parabolic  differential  equations. 
Similar  comments  apply  to  problem  (6.3)  -  (6.6)  when  f(x)  vanishes  and 

the  boundary  values  gQ  (t)  or  g1  (t)  are  non- zero.  We  are  thus  led 

to  consider  the  following  two  special  forms  of  problem  (6.3)  -  (6.6) 


and 


h 

-  aD2| 

|  u  =  0 

0  < 

t 

<  T, 

u 

(x,Q) 

=  h 

n 

(x-x  )  0 

o' 

< 

X,  3 

u 

(0,t) 

=  0 

0 

< 

t  _< 

u 

(*.t) 

=  0 

0 

< 

t  £ 

[D< 

-  aD2 
:  xj 

|  w  =  0 

0 

< 

t  £ 

w. 

(x,0) 

=  0 

0 

< 

X  £ 

w 

(0,t) 

=  h 

Y 

(t)  •  0 

< 

t 

w 

(Z,  t) 

=  0 

0 

< 

t  £ 

o  — 


(BVP). 


(BVP) 


In  (BVP)  we  require  that  n  is  a  non-negative  interger;  in  (BVP)  v 
o  4 

can  be  any  non-negative  real  number.  A  solution  to  these  problems  can 
be  obtained  by  defining  an  extension  of  the  initial  values  to  the  entire 
real  line  in  such  a  way  that  the  solution  of  the  resulting  Cauchy 
problem  must  also  satisfy  the  correct  boundary  conditions  along  x  =  0 
and  x  =  i 

For  problem  (BVP)  consider  the  following  periodic  extension  of  the 
given  initial  values : 
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£  (x) 


hn  (x  -  2k£  -  x  )•,  2kS,  <  x  <  (2k  +  1)  Si 

hn  (x  -  (2k  +  2)1  +  xQ),(2k+l)&  lx  <  (2k  +  2)1 


where  k  =  0,  ±1,  ±2,  ...  .  This  function  repeats  with  a  period  2£  and 
is  antisymmetric  about  any^  of  the  points  x  =  k&,  where  k  =  0,  ±1,  ±2, 
...  .  A  single  £ycle  of  f  is  indicated  in  Figure  1  over  the  interval 
[0,  2Jt] .  Since  f  is  antisymmetric  about  x  =  0  and  x  =  l  in  particular 
then  the  solution  of  the  Cauchy  problem 


D-D 
t  x 


u  =  0 

u(x,  0)  =  f  (x) 


(IVP) 


will  also  be  antisymmetric  about  x  =  0  and  x  =  &  and  must  therefore 
be  a  solution  to  (BVP) ^  as  well.  To  obtain  the  solution  to  (IVP)3  in 

a  useful  form  note  first,  that  in  the  interval  [0,  2 £]  the  function  f 
can  also  be  written  as  the  sum  of  jump  functions  in  the  form 

n 

f  (x)  =  hn  (x-xo)  +  ^  £k  hk  Cx'*°  +  C"1)Ilhn  (X'2*  +  V>  C6-7;) 

k=o 


o  <_  x  2  where  f  denotes  the  jump  in  the  k-th  derivative  values 

K. 

of  f  at  x  =  H,  namely  , 

_  x  )n_k 

fk  '  'Cl  +  C-D*)  (n-  kjl  •  C6-« 

By  adding  on  all  other  jumps  occuring  to  the  right  and  left  of  [0,2£] 
in  proper  sequence  we  obtain 

oo  n 

£  (x)  =  ^T^jhn  (x-2j &-xq)  fk  hk  ((x-(2j+l)A) 

j'=o  k=o 

+  C-l)n  hn  (x-  (2j  +2) it  +  XQ)J  (6.9) 

°°  n 

Z  |hn  *  V*  L  fk  h*  (X+  (2j  +1)A) 

j=o  k=o 


+  (-l)n  h*  (X  + 


( 2 j  *  2)1  -  xo) 


] 
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Figure  1.  Graph  of  f(x)  in  the  interval  [0,2fc]. 
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The  solution  to  (IVP)^  can  thus  be  written  (formally  at  least) 


+  (-1)  Hn  (x- (2j+2) £+xr 


(x,t)  =^[^Hn  (x-2jA-xQ,at)  +^fk  Hk  (x-(2j+l)A,at) 
i=o  k=o 

o'at)] 

H*  (x+2jA+xo,at)  +  fR  H*  (x+(2j  +  l)A,at) 

=  0 

-X0’at)J 


qo 

St 


+  (“I)  H*  (x+(2j+2) i- 


(6.10) 


We  denote  the  truncation  of  this  to  N  terms  by 

N  n 

Vx>t)  =S  Hn  C^-2jit-xo,a.t)  fk  Hk  (x-  (2j +1)  SL, at) 

j  =o  k=o 

+  C“l)n  Hn  (x- (2 j +2) &+xq , at)  1 

N  n  -* 

-  E[Hn*  (x+2jA+xo,at)  +2fkHk  Cx+(2j+l) A,at)J 


+  C-l)n  H*  (x+(2j  +  2)  £-XQ,at) 


(6.11) 


The  accuracy  of  u  as  an  approximation  for  u  in  the  domain  [0,£]  x 
[0,T]  can  be  estimated  by  comparing  the  values  of  u^.  along  x=0,  t=0 

and  x=£  with  the  prescribed  values  for  u.  We  have 


h  (x-x  )  n  t=0 

n  °  t — ' 

H*  ((2N+l)A+xQ,at)  +  ^fR  H*  ((2N+2) A,at) 

k=o 

+  (-l)n  H*  ((2N+3)A-xo,at),  x=£ 

It  is  clear  that  u^  satisfies  the  boundary  and  initial  data  exactly 

along  x=0  and  t=0  but  does  not  vanish  identically  along  x-H,  as 
required.  However  using  (4.16)  the  truncated  series  can  be  bounded 
along  x=i  in  the  form 

iun t)  |  1  const,  exp  (xq+  (2N+1) JL) 2/4at  j  . 


uN(x,t)  = 
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Since  u  and  are  both  solutions  of  (6.3)  then  the  maximum  principle 
for  parabolic  differential  equations9  states  that  their  maximum  difference 
occurs  on  the  boundary.  Thus  the  solution  u  of  (BVP)  can  be  expressed 
in  the  form  * 


u(x,t)  =  u^x^)  +  (x,t) 


(6.12) 


where  is  a  remainder  term  satisfying 

|^(x,t)|  <_  const,  exp  j^-(xo+(2N+l)5,)2/4atJ 

uniformly  in  [o,£]  x  [o,T].  The  approximation  u,  for  u  is  therefore 
particularly  accurate  for  small  values  of  t.  In  fact  for  many 
applications  it  is  sufficient  to  use  only  the  lowest  order  approximation 
uq  (x,t). 

A  similar  analysis  applys  to  problem  (BVP)  which  has  the  equivalent 
initial  value  problem  4 

-  -  aD2l  w 

t  X 


[ut  -  <1 

w(x,0)  =  2a'1  y~||h«;v  fx*2]2j  -  h2,f(x-(2ja2jijl 


(I  VP) 


with  the  formal  solution 


w(x,t)  =  2a“Y  ^  ^ 


In  this  form  it  is  apparent  that  w(x,t)  is  antisymmetric  and"* thus 
vanishes  across  x=l.  On  the  other  hand,  by  rearranging  the  terms  into 
the  equivalent  form 

-Y 


H2Y  (x+2j£,at)  -  (x-  (2j  +  2)  £,at) 


(6.13) 


w(x,t)  =  2a"'  H*Y(x,at) 


+  2a 


•yE 

j=l 


H2Y(x+2j&,at)  -  H2y  (x-2j l, at) 


(6.14) 


it  becomes  clear  that  w(x,t)  satisfies  the  prescribed  boundary  condition 


along  x=0.  Denoting  the^truncated  series  by  ' 


wN(x,t)  =  2a"Y^H*Y(X+2j£,at)  -  H^x-  (2j+2)  4,at)  (6.15) 

u„.  _j=°4 


we  can  show  as  in  the  previous  case  that 


w(x,t)  -  wN(x,t)  +  R^x.t)  (6.16) 

9/1  *  Tikhonov  and  A.  A.  Samarskii ,  "Equations  of  Mathematical 

Physics" ,  Page  206,  Pergamon  Press ,  Inc.,  MacMillan  Company, 

New  York,  1963. 
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where  _  r 

|l^(x,t)|  <_  const  exp  I- ((N+2)  £)  "/4at 

uniformly  in  [0,£]  x  [0,T]. 

To  conclude  this  section  let  us  briefly  touch  on  two  problems  where 
the  function  derivative  (flux)  values  are  specified  at  the  boundary. 
Consider 


[Dt  -  aDx 

u(x,o)  =  hY(x-XQ) 
ux(0,t)  =  0 
ux(*,t)  =  0 

and 

jl)t  -  aD^  j  w  =  0  (BVP)  6 

w(x,o)  =  0 
wx(0,t)  =  h  (t) 

Wx(£,t)  =0 


By  derivations  similar  to  that  for  (BVP)^  and  (BVP)^,  only  using 

symmetry  in  place  of  antisymmetry,  we  can  arrive  at  the  following 
solutions 

™  n 


u(x,t)  = 


Ek 

j=oL 

Ek 


(x-2j£-xQ,at)  -  fR  Hk(x-(2j+l)£,at) 
k=o 


(6.17) 


("l)n  H  (x-(2j+2)£+x  ,at) 


E 


(x+2jJl+xo,at)  -  fk  H*k  (x+(2j+l)A,at) 
k=o 


(-l)n  H*  (x+(2j+2)£-xo,at) 


and 


w(x,t)  =  2a  ^ 
for  (BVP)^  and  (BVP)^  respectively. 


'£  h- 


(x+2ji,at)  +  H2y+1  (X- (2j+2) Z 


,at)J  (6. 


18) 
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We  should  mention  that  the  conversion  of  initial-boundary  value 
problem  to  equivalent  Cauchy  problems  works  only  for  the  diffusion 
equation  and  not  for  more  complicated  parabolic  equations  without 
special  modification.  Nevertheless  in  these  cases  we  can  use  the 
Diffusion  Equation  Solution  Sequence  method,  alluded  to  above,  to 
obtain  approximations  which  are  comparable  with  the  first  term  of  the 
expansions  just  derived.  This  technique  has  been  discussed  in  detail 
for  Dirichlet  type  of  boundary  conditions  (function  values  specified) 
and  in  a  later  report  will  be  applied  to  Robin's  type  or  connective  heat 
transfer  boundary  conditions,  such  as  occur  in  gun  barrels. 


VII.  CONCLUSION 


The  function  H^,  H*  and  v^  have  been  defined  and  their  properties 

investigated.  In  particular  they  have  been  shown  to  be  solutions  of 
the  diffusion  equation  with  the  special  initial  values 


where 


and 


H 

Y 

(x,0) 

II 

-< 

/ — s 

H* 

Y 

(x,Q) 

=  h*  (x) 

V 

n 

(x,0) 

=  xn/n! 

(  0  x<0 

h 

(x)  = 

Y 

{  xr/y\  x>0 

h* 

Y 

(x)  = 

\(-x) 

The  functions  and  H*  were  then  used  to  develop  the  series 

expansions  (6.10),  (6.13),  (6.17)  and  (6.18)  for  the  solutions  of 
initial-boundary  value  problems  (BVP)^  -  (BVP)^. 


\ 
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APPENDIX  THE  POLYNOMIALS  p  AND  q 

o,n  no,n 

n  P«  r,  Cx)  =  v  fx) 

_  o,n  v  _ n  _ 

<o  o 

0  1 

1  x 

2  (x2+2)/2 

3  (x3+6x)/3! 

4  (x4+12x2+12)/4! 

5  (x5+20x3+60x)/5: 

6  (x6+30x4+180x2+120)/6.' 

7  (x7+42x5+420x3+3360x2+840x)  /7 

8  (x8+56x6+840x4+3360x2+1680)/8l 

9  (x9+72x7+1512x5+10,080x3+15,120x)/9.' 

10  (x10+90x8+2520x6+25,200x4+75, 600x^+30, 240) /10J 

11  (x1;L+110x9+3960x7+55,440x5+277,200x3+332,640x)/11] 

12  (x12+132x10+5940x8+110,880x6+831,600x4+l,995,840x2+665,280)/12i 


In  general 


P 

o,n 


(x) 


V 

n 


(x)  = 


0  n  =  -1,-2 , . . . 

[n/2] 

Exn"2k 

(n-2k) Ik  I  n  =  0,1,2,... 
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n 


q  (x) 
o,n  v  J 


0  0 

1  2 

2  x 

3  2(x2+4)/3! 

4  ,  2(x3+10x)/4! 

5  2(x4+18x2+32)/5! 

6  2(x5+28x3+132x)/6! 

7  2(x6+40x4+348x2+384)/7 

8  2 (x7+54x5+740x3+2232x) / 8 1 

9 .  2  (x8+70x6+1380x4+7800x2+5568)  /9 

10  2  (x9+88x7+2352x5+21 , 120x3+45, 744x)/10.' 

11  2 (x10+108x8+3752x6+48, 720x4+201, 744x2+111 , 360) /II ! 

12  2  (x11+130x9+5688x7+100,464x5+666,  384x3+1, 117, 728x)/12.' 


In  general 


[Cn-l)/2] 


q  (x) 
no,nv  J 


where  for  n=l,2,3, 


nk 


=  (2/n!)  ^ 


a  ,  x 
nk 


n-l-2k 


k=o 


1  for  k=o 

Vl,k  +  2(n‘1)  an-2,k-l  £or  k=1>2‘ 


[01)72] 
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n  q  (x) 

 o,n   ' 

-1  1 

-2  -x/2 

-3  (x2-2)/ (-2) 2 

-4  (x3-6x)/C-2)3 

-5  (x4-12x2+12) / (-2) 4 

-6  (x5-20x3+60x)/(-2)5 

-7  (x6-30x4+180x2-120)/(-2)6 

-8  (x7-42x5+420x3-840x)/(-2)7 

-9  (x8-56x6+840x4-3360x2+1680)/ C-2)  8 

-10  (x9-72x7+1512x5-10,080x3+15, 120x)/ (-2) 9 

-11  (x10-90x8+2520x6-25,200x4+75,600x2-30,240)/(-2)10 

-12  (x11-110x9+3960x7-55,440x5+277,200x3-332,640x)/(-2)11 


In  general 


(-2)1_n 


[  (n-l)/2] 

E 


k=o 


(-Dk(n-1)! 

(n-l-2k) I  k! 


l-2k 
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